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Abstract 

Let n C 1" be a C 2 bounded domain and x > 0, A > 0, be constants. We 
obtain upper bounds for the quenching time of the solutions of the nonlocal 
parabolic MEMS equation u t = An + A/(l - u) 2 (l + X J n 1/(1 - u)dx) 2 in 
x (0, oo), u = on dQ, x (0, oo), u(x,Q) = uq in C2, when A is large. We 
will prove the compactness of the quenching set under a mild condition on 
uq. When f2 = Br = {x £ R" : \x\ < R} and uq is radially symmetric and 
monotone decreasing in < r < R, we prove that the point x = is the only 
possible quenching set. When f2 = Br and uq is radially symmetric which also 
satisfies some strict concavity assumption, we prove that for any f3 € (2, 3) the 

solution satisfies 1 — u(x,t) > C\x\p for some constant C > and that the 
solution u quenches in a finite time for any sufficiently large A > 0. We also 
obtain the quenching time estimate in this case. 
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Micro-electromechanical systems (MEMS) devices are key components of many 
electronic devices including accelerometers for airbag deployment in cars, inkjet printer 
heads, and the device for the protection of hard disk, etc. It is therefore interesting to 
understand the mathematical modelling of the MEMS devices and study the various 
properties of such models. Interested readers can read the book, Modeling MEMS 
and NEMS [22J, by J. A. Pelesko and D.H. Berstein for the mathematical modeling 
and various applications of MEMS devices. 
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One model of MEMS |10| . [19] . [22] . consists of a fixed membrane and a deformable 
membrane of the same shape which is coated with a thin dielectric material and 
held fixed at the boundary. When no voltage is applied to the membranes, the two 
membranes are parallel to each other with one at a fixed small distance on top of 
the other. When a voltage is applied to the membranes, the deformable membrane 
will deflect towards the fixed membrane. Recently there are a lot of study on the 
equations arising from such model of MEMS by N. Ghoussoub, Y. Guo, Z. Pan and 
M.J. Ward [3], 0], [5], ©, [9], [10], J.S. Guo, B. Hu and C.J. Wang [8], K.M. Hui [H], 
[15], p], N.I. Kavallaris, T. Miyasita and T. Suzuki [TF], F. Lin and Y. Yang [19], 
L. Ma, Z. Guo and J.C. Wei p], P2], [20], G. Flores, G.A. Mercado, J.A. Pelesko 
and A. A. Triolo [2], [21], [23] etc. 

As observed by N. Ghoussoub, Y. Guo, J.A. Pelesko and others [13], [22J, the 
deflection of the deformable membrane from its equilibrium position is modeled by 
the following parabolic equation, 



du 
dt = 
u =0 
k u(x,0) =u 



Au + 



A/(x) 
1-uf 



in Q x (0,T) 

on dQ x (0,T) 
in Q 



where A > is a constant and / ^ is a nonnegative function on Q which depends on 
the dielectric constant of the coating on the membrane. When the voltage between 
the membranes are due to circuit series capacitance, the deflection of the deformable 
membrane from its equilibrium position is modeled by the following nonlocal parabolic 
equation [TJ, 



du 
~dt 



--Au + 



A 



:i--) 2 (i+xLi^y) 2 



u =0 

k u(x, 0) =Uq 



in tt x (0,T) 

on dQ x (0,T) 
in Q 



where A > and x > are a constant. In [19] F.H. Lin and Y. Yang by usin^ 
variational argument derived the following nonlocal MEMS equation 



-Av 



A 



in Q 



v =0 on dQ 



for modeling the stationary deflection between the two parallel plates of an electro- 
static MEMS device with circuit series capacitance. 

Note that A is proportional to the square of the voltage. When the voltage is very 
large, the two membranes will touch each other or quench and we would expect the 
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solutions of (PQ and (P\) cease to exist after a finite time. One would like to get 
estimates for the touchdown time of the solutions of (Pa) and (P{). We refer the 
readers to the papers [4], [5], [6], by N. Ghoussoub and Y. Guo for various touchdown 
time estimates for the solutions of (P'x)- 

In this paper we will obtain touchdown time estimates for the solutions of (-Pa)- 
We prove that the quenching set of (Pa) is compact under a mild condition on the 
initial value of the solution. When Q = Br = {x G M n : \x\ < R} and Uq is radially 
symmetric and monotone decreasing in < r < R, we prove that the solution u of 
(Px) satisfies 

1 -u(x,t) > C\x\ 2 in B R x[t ,T) (1) 

for some constants t G (0,T), C > 0, and the point x = is the only possible 
quenching set. When Q = Br and Uq is radially symmetric which also satisfies some 
strict concavity assumption, we prove that for any (3 G (2, 3) the solution u of (Pa) 
satisfies 

1 -u{x,t) > C\x\? inQl (2) 

for some constant C > and that the solution u quenches in a finite time for any 
sufficiently large A > 0. We also obtain the quenching time estimate in this case. 

We will assume that Q C M ra is a bounded C 2 domain and A > 0, x > 0, for the 
rest of the paper. We start with some definitions. For any 5 > 0, R > 0, T > 0, let 
fi^ = {x G : dist(x, dft) < 5} and = Br x (0, T). Observe that when Q is a 
bounded convex domain then there exists a constant Sq > such that for any x G Qs 
there exists a unique point y G dfl such that the line seqment xy is perpendicular to 
dfl at y. 

For any constants x > 0, A > 0, and 

u G P 1 (fi) with Uq < a a.e. in Q (3) 

for some constant < a < 1 we say that u is a solution (subsolution, supersolution 
respectively) of (P A ) in Vt x (0,T) if u G C 2 '\n x (0,T)) n C(Q x (0,T)), u < 1, 
satisfies 

du A ^ . m , 

— = Au + in ft x 0, T) 

dt (i--) 2 (i + xLt^) 2 

(<, > respectively) in the classical sense with u(x,t) = (<, > respectively) on 
dn x (0,T), 

sup u(x,t) < 1 V0 < T < T, 

Qx[0,T>] 

and 

\\u(-,t) - MqIU 1 ^) ~* as t 0. 
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For any solution u of (Pa) we define the quenching time or touchdown time T\ > 
as the time which satisfies 

supu(x, t) < 1 VO < t < T x 
n 

lim sup u(x, t) = 1. 

We say that u has a finite quenching time if T\ < oo and we say that u quenches at 
time infinity if Tx = oo. For any solution u of (P\) we let the quenching set of u to 
be the set of points iGfi such that there exists a sequence (xk, t k ) G Q x (0, T A ) such 
that Xfc — )■ x and u(xk, tk) — > 1 as k — > oo. 
We first recall some results of [TH] and [To] . 

Theorem 1. (Theorem 2.1 and Theorem 2.2 of \ffl$) Let 

— b < u Q < a in Q (4) 

for some constants < a < 1 and 6 > 0. T/ien for any A > and x > t/iere exists 
T > suc/i i/iai (Ta^ a unique solution —b<u<linQx(0,T) which satisfies 

u{x,t)= [ G(x,y,t)u (y)dy + X [ [ ^77^ — ^ Tz — rzdyds (5) 

Jn Jo Jn (1 - u(y, s)) 2 (l + X J n j^F) ) 

for all (x,t) G Q x (0,T) where G(x,y,t) is the Dirichlet Green function for the heat 
equation infix (0, T). 

Lemma 2. (c/. Theorem 2.1 of \W$ ) Let ^0,1,^0,2 £ 6e such that < «o,i < 

W0.2 < a in O /or some constant < a < 1. Let A > and < / £ C(fi x (0, T)) n 
L°°(f2 x (0, T)). Suppose u\, ui, are nonnegative subsolution and supersolution of 
(P' x ) in Q x (0, T) with initial value Uq = Uo,i,Uo,2, respectively. Then u\ < U2 in 
Ox (0,T). 

By an argument similar to the proof of Proposition 2.1 of [5] we have the following 
theorem. 

Lemma 3. Let Vt d W 1 be a bounded convex domain. Let uq satisfy (4) for some 
constants < a < 1 and b > 0. Suppose there exists a constant 5q > such that 

u o\n 6 e C 1 (^<5o) an d ~q^~ < on 90 (6) 

where djdv is differentiation with respect to the unit outward normal v on dfl and u 
is the unique solution of (P\) in Q x (0, T) given by Theorem 1. Then there exists a 
constant S\ G (0, So) such that for any x G if y is the unique point on dQ such 
that the line segment xy is perpendicular to dfl at y, then 

d 

—u{z,t)<0 Vzexy,Q<t<T 
where d/dn is differentiation with respect to the unit vector n along the direction xy. 
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By Lemma 3 and an argument similar to the proof of Proposition 1.3 of |T6] and 
Propsotion 2.1 of [9] we have the following two corollaries. 

Corollary 4. Let Q C R n be a bounded convex domain. Let uq satisfy (4) and (6) 
for some constants < a < 1, b > 0, and 5 > 0. Suppose u is the unique solution 
of (P\) in Q x (0,T) given by Theorem 1. Let 5\ be as in Theorem 3. Then there 
exist a\ > and 5 G (0, 5\) such that for any y £ Q$ there exists a fixed-sized cone 
T(y) C with vertex at y such that \T(y) \ Qs\ > ai and u(z,t) > u(y,t) for any 
z GT(j/) and < t < T. Moreover 



Corollary 5. Let Q C W 1 be a bounded convex domain. Let uq satisfy (4) and (6) 
for some constants < a < 1, b > 0, and Sq > 0. Suppose u is the unique solution of 
(P\) in Q x (0, T) given by Theorem 1 such that u touchdown at time T. Then the 
set of touchdown points for u is a compact subset of Q. 

By (7) and an argument similar to the proof of Theorem 4.3 of [16J we have the 
following theorem. 

Theorem 6. Let Q C R™ be a bounded convex domain. Let uq satisfy (4) and (6) 
for some constants < a < 1, b = 0, and Sq > 0. Let \ > 0. Then there exists a 
constant C\ > and such that for any X > Xq = C\yu\ and any solution u of (P\), u 
quenches in a finite time 



Lemma 7. Let X > 0, x > 0, and f2 C M™ be a bounded domain such that Br C Q 
for some constant R > 0. Let < u < a in Q for some constant < a < 1 and let 
< u < 1 be a solution of (P\) in Q x (0, T) with 




V0 < t < T. 



(7) 




Let 





and < S\ < inf <t<T A(i). Suppose X > 2n/5iR 2 . Then T satisfies 



Proof: Let v = 1 — u and vq = 1 — Uq. Then v satisfies 

v t = Av- \A(t)v~ 2 in tt x (0, T) 
f(x,t) = l on<9fix(0,T) (10) 

v (x, 0) = Vq in Q. 

Let 



where 



Then 



ip(x,t) = 1 — X5ic tyl 
2n 

c = " 



R 2 



A^i? 2 ' 

2n 1 



# 2 A(5ic 



en 1 - 1. (11} 



Hence by (11), 



ip t - Aip + \5xip~ 2 > A5ic [-1 - (2nt/R 2 ) + c l ] > in 
where T\ = min(T, l/(A<5iC )). Hence ip is a supersolution of 

ipt = Aip - \5 1 tp~ 2 in Qg 
V> = 1 on x [0, Ti) nB R x {0}. 

Suppose T > l/(A<5iC ). Since v is a subsolution of (12), by Lemma 2, 



(12) 



< v < ip in 
=> < v{0, l/(A5 lCo )) < ^(0, l/(A5 lCo )) = 0. 

Contradiction arises and (9) follows. □ 

Theorem 8. Let < Uq < a be a radially symmetric function in Br C W 1 for some 
constant < a < 1 which is monotone decreasing in < r < R. Let T > 0. Suppose 
u is the unique solution of (P\) in given by Theorem 1 and u touches down at 
time T. Then there exist constants to G (0,T) and C > depending on A and x such 
that (1) holds. Hence x = is the only quenching point of u at time T . Moreover for 
n > 3 we have 

sup / i — < 00 • ( 13 ) 

o<t<Tj BR 1 - w(M) 

Proof. We first observe that by an argument similar to the proof of Theorem 4.3 of 
[To] and Theorem 1.5 of [13] 1 > u(x,t) = u(\x\,t) > is radially symmetric in 
with 

M r (0, t) = and u r (r, t) < V0 < r < i?, < t < T. (14) 
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Let 

Vl (x,t) = I G(x,y,t)u (y)dy 



f B R 

and 

F(x, t,s) = X I G(x, y, t - s)A{s){l - u(y, s))~ 2 dy 

JB R 

where A(t) is given by (8) with Q = Br. Since < v 1 < a satisfies the heat equation 
in with v% = on OBr x (0,T) and u$ is radially symmetric and monotone 
decreasing for < r < R, by the maximum principle and an argument similar to 
the proof of Theorem 1.5 of [13] V\ is radially symmetric and t>i jr (r, t) < for any 
< r < R and < t < T. 

Similarly for any < s < T, F(-, •, s) is a radially symmetric solution of the heat 
equation in Br x (s, T) with F(x, t, s) = for any x G <9I?r, < s < t < T, and 

OF 

— (r,t,s)<0 VO < r = d < i?, < s < t < T. (15) 
or 

By Theorem 1 u satisfies (5). By (5) and (15), 

u r (r,t) =v lir (r,t) + / F r (r,t,s)ds< j F r (r,t,s)ds VO < r < R, T/3 < t < T. 
Jo Jo 

(16) 

Let 

T 

v 2 (x,t) = 3 F(x,t,s)ds V|x| < R,T/3<t <T. (17) 
Jo 

Then v 2 is radially symmetric and satisfies 

v 2 ,t = Av 2 = -^(r™- 1 ^,,), in B R x (T/3, T) 

w 2 = on dB R x (T/3, T) ( 18 ) 

w 2 > iiiB B x {T/3} 

By (15) for each T/3 <t <T V2 t r{ r ,t) is monotone decreasing in < r < R. Hence 
for each T/3 <t <T v 2 (x,t) will attain its maximum at x — 0. Differentiating (18) 
with respect to r and letting q = r n ^ 1 v 2ir , 

Ti — 1 

Qt = Qrr Qr ™ (0, R) x (T/3, T). (19) 

r 

Let q = q + er n for some constant e > to be determined later. Then 

Ti — 1 

qt - q rr + q r = in (0, R) x (T/3, T) (20) 

r 

and 

g(0,t)=0 VT/3<t<T. (21) 
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By the maximum principle ([T], [18]) there exists a constant C\ > such that 

v 2 , r (R, *) < ~Ci VT / 2 <t<T. 

Hence by choosing < e < e\ = C\j (2R) we have 

q{R, t) < -(C 1 /2)R n ~ 1 VT/2 < t < T. (22) 

By the strong maximum principle, 

v 2 (r, t') < max v 2 {y, t) = v 2 (0, t) VO < r < R, T/2 < t < t' < T 

v 2 (0, t') < v 2 (0, VT/2 < t < t' < T (23) 

f2,t(0,t)<0 VT/2<t<T. 

Suppose i?2,t(0, t) = on (T/2, T). Then 

u 2 (0, t') = u 2 (0, t) VT/2 < t < t' < T. 

This contradicts (23). Hence there exists to £ (T/2,T) such that 

MMo)<0. (24) 

By (18), 

?2 — 1 

V2,t(r, t ) = ^2,rr(?", *o) H «2,r(»", to)- (25) 

r 

Since t>2, r (0,t) = for any T/3 < t < T, letting r — > in (25) by the l'hosiptal rule, 

U2,i(0, to) = nu 2 ,rr(0, t ). (26) 

By (24) and (26), 

f2,rr(0,t ) < 0. 

Then there exists < r% < R such that 

e 2 = — max v 2>rr {r, t ) > 0. (27) 

0<r<ri ' 

Let £3 = — max ri < r <fl v 2tr (r, to). Then e 3 > 0. Let £0 = min(ei, £2/2, e^/{2R)) and 
< e < £q. Since t>2,r(0, to) = 0, by the mean value theorem for any < r < r\ there 
exists < r' < r such that 

q(r, t ) = r"( ^ r( ^ to) + e) = r n (v 2 , rr (r', t ) + e) < 0. (28) 

and 

g(r, t ) = r n " 1 (v 2 ,r(r', t ) + er) < r n_1 (-£ 3 + eR) < Vn < r < i2. (29) 
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By (28) and (29), 

g(r,t )<0 V0<r<#. (30) 
By (20), (21), (22), (30) and the strong maximum principle, 

q(r, t) < V0 < r < R, t < t < T 

V2, r (r, t)<-er V0 < r < i?, t < t < T (31) 

Hence by (16) and (31), 

u r (r, t) < v 2 , r (r, t) <-er V0 < r < R, t Q < t < T. 

Thus 

pr pr 

u(r,t) = u(0,t) + / Ur(p,t)dp<l-e pdp < 1 - (e/2)r 2 
Jo Jo 

holds for any < r < R, to < t < T, and (1) follows. By (1) x = is the only 
quenching point of u and for any n > 3, 



< C / r n - A dr < CR n - 2 Vt < s < T 



Br 1-U{Z,S) jo 

and (13) follows. □ 

Theorem 9. Let < u < a be a radially symmetric C 2 function on Br C M. n for 
some constant < a < 1 which is monotone decreasing in < r < R with 

u = on dB R and <(r) < -ci V0 < r < i? (32) 

/or some constant c\ > 0. Lei x > and T > 0. T/ien there exists a constant A > 
such that for any A > A and 2 < f3 < 3 if u is the unique solution of (P\) in 
given by Theorem 1, then there exists a constant C > such that (2) holds and for 
any n G Z + (13) holds. Moreover there exists a constant \\ > Xq such that if X > Ai 
and T is the maximal time of existence of the solution u of (P\), then T satisfies (9) 
and x = is the only quenching point of u at the quenching time T . 

Proof. Note that the estimate (2) for the case n = 1 and sufficiently large A is proved 
in [5]. Here we will modify the proof of [8] so that it works for any dimension n. Let 
A > A and 2 < (3 < 3 for some constant A > to be determined later. Let u be the 
unique solution of (P\) in given by Theorem 1. By an argument similar to the 
proof of Theorem 8 < u < 1 is radially symmetric in and u satisfies (14). Let 
v — 1 — u and v o = 1 — «o- Then v satisfies (10). By (32), 

v = 1 on dB R and v' (r) > c x r V0 < r < R. (33) 
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By (10), (33), and an argument similar to the proof of Proposition 2.1 of [9 J and 
Theorem 4.1 of [8] there exists a constant < c 2 < C\ depending on u but is 
independent of u such that 

v r (R,t)>c 2 R V0<t<T. (34) 
Let w = v 13 . Then w satisfies 

w t - ^=- x {r n - l w r ) r = -\^A{t)w^ - ^ • ^ in (0,R) x (0,T) (35) 

with «j = lon dB R x (0, T). By (14) w r (r, t) > for any < r < R and < t < T. 
Differentiating (35) with respect to r and setting q = r n w r , we have 

&-g rr + g r =A(3 - /3)A(t)tu ?g + ^— r n ^-2-^ 



r p p w 

>A(3 - g - 2 ^ ~ ^ —q in (0, R) x (0, T). (36) 

p w 



qyg n — 1 2 



Now 

W «-<? = (^l) 1=^1- '—^^ ^ QrWr- (37) 

Since < w < 1, by (36) and (37), 



ft - ftr + —ft >- f A(3 - /3)A(t)r»" V-f - 2^-^g f 
r u> \ p 

>^ (a(3 - /3)A(t)r"- 1 - 2^-^ft) in (0, it!) x (0, T). (38) 

We now choose e > sufficiently small such that 

A(0) > S 1 := [1 + X |w n _ 1 |(2/e)7(7i - (2/ f3))- 1 R n ^]- 2 (39) 

and 

e < /3min(ci(l -o)^ -1 ,ca) (40) 
where is the surface area of an unit ball in IR n . We next choose 

2en(/?-l) 

An > 



/3(3 - /3)5 1 

Since A > A , 



A(3 - P)8 X > 2 ^ l \ n. (41) 
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By (39) and (41) there exists < 7\ < T such that [0, Ti) is the maximal interval 
such that 

A(3 - 0)A(t) > 2 ^~ 1 h n (42) 

holds for any < t < T x . Let q = q - er n . Then by (38) and (42), 

n — 1 (6 — 1) w 

qt ~ q rr + q r + 2 V ' \ J — q r > in (0, R) x (0, 7\). (43) 

r p w 

Now by (34) and (40), 

q(R,t) = R^ipv^VriR^) - sR) = R^ipVriR^) -eR)>0 V0 < t < T. (44) 
By (33) and (40), 

q(r, 0) > V0 < r < R. (45) 
By (43), (44), (45) and the maximum principle, 

q(r, t) > V0 < r < i?, < t < T x 
=> w r {r, t)>er V0 < r < R, < t < Ti 

=>- w(r, i) > w(0, t) + E -r 2 > £ -r 2 V0 < r < R,0 <t < T x 

1_ 

^ f(r,t)>(-) A V0 < r < R,0 < t < 7\. (46) 



By (46), 



/7 y i o 

< |o; n _ 1 |(2/e)?(n-(2//3))- 1 J R n ^ V0 < t < ^ (47) 



v(z,t) 

A(t) > 5i V0 < t < Tx. (48) 

If T\ < T, then by (41) and (48) there exists T 2 G (Ti,T) such that (42) holds on 
(0,T 2 ). This contradicts the maximality of T ± . Hence Ti = T. By (46) and (47) we 
get that (2) and (13) hold. 

We now let A > Ai for some constant Ai to be determined later and let T be the 
maximal time of existence of the solution u of (-Pa)- Suppose T > 1. By repeating 
the above argument but with T being replaced by 1 in the argument and noting that 
then the constant c 2 is now independent of T we can find constants e, 5i, and Ao as 
before which are then all independent of T. Then (2) holds in and (13), (48), 
hold for any < t < 1. Let 

Ai = max(A , An/^R 2 ), 3/<?i). (49) 

By (49) and Lemma 7, 



Contradiction arises. Hence T < 1. Note that then we can choose the constant C2 in 
the above argument to be independent of T. Hence the constants e, Si, and A are 
all independent of T. Then by (48) Lemma 7 holds. Hence T satisfies (9). Thus u 
quenches in a finite time T bounded above by the right hand side of (9). By (2) x — 
is the only quenching point of u at the quenching time T and the theorem follows. □ 
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